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(~| ■ Abstract 

Infinitely many particles of two types ("plus" and "minus") jump randomly 
along the one-dimensional lattice = eZ. Annihillations occur when two particles 
of different time occupy the same site. Assuming that at time f = all "minus" 
particles are placed on the left of the origin and all "plus" particles are on the right 
of it, we study evolution of /3s{t), the boundary between two types. We prove that 
^ in large density limit e — )• the boundary /^^(t) converges to a deterministic limit. 

, This particle system can be interpreted as a microscopic model of price formation 

' on economic markets with large number of players. 

^ ! Keywords: stochastic particle systems with annihilation, scaling limits, microscopic 

^ I models of price formation 

(N ■ MSG classes: 60J99, 60K35, 91B26 
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1 Introduction 

On one-dimensional lattice = eZ = {em : m e Z} , e > 0, there are particles of 
two types — "plus particles" and "minus particles". Denote by i^^it) the number of 
plus(minus)-particles at site em at time t. We define a continuous time Markov process 
on [0, oo) by the following conditions: 

(1) at time all plus particles have positive coordinates, all minus particles have nega- 
tive coordinates; 
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(2a) any plus particle, independently of other particles, performs a simple random walk: 
that is it jumps from em to e{m + 1) with rate /i+ and from em to e{m — 1) with 
rate A+; 

(2b) any minus particle, independently of other particles, performs a simple random 
walk: that is it jumps from em to e{m + 1) with rate A_ and from em to e{m — 1) 
with rate 

(3a) if a plus particle jumps to a site where there are minus particles it immediately 
annihilates with one of the minus particles at this site; 

(3b) if a minus particle jumps to a site where there are plus particles it immediately 
annihilates with one of the plus particles at this site. 

At any time t > the state of the process is the vector (z/^(t),m G Z). However, it 
follows from (3a) and (3b) that for any m and t 

Moreover, all minus particles are always to the left of the leftmost plus particle. It will 
be convenient to define /3e(t) G as the point where the last annihilation before time t 
happened. We call it the phase boundary. 

Note that there are no problems with the existence of this process. 

Besides the interpretation related to annihilation of particles there is another one — 
the microdynamics of the price formation, where the market contains many players and 
is formed by their behaviour. Namely, /3e(t) is the price of some product at time t. 
Plus particles (bears) want to lower the price of this product (we assume further on 
that = A+ — > 0), minus particles (bulls) want to increase the price (we assume 
a_ = A_ — /i_ > 0). Annihilation is a bargain which is performed when the demand 
and offer prices meet together. Recent models of price formation [5], El [3l H] have much 
in common with our model, however they are closer in spirit to queueing models. Our 
model is closer to statistical physics models. Anyway, all such models cannot pretend on 
practical implementation, mainly because external influence on the action of players is 
not taken into account. 

We consider the large density limit e — )■ under the time scaling t = Te~^, where t 
is microtime and r is macrotime. Our goal is to find asymptotic behaviour of the price 
/3(r) = lims^o f3e(t) as the result of many micro-bargains. 

2 Main result 

Initial distribution of particles 

We assume for simplicity that at time t = the distribution of plus particles is (inho- 
mogeneous) Poisson with density p+{em), where p+{x) is some strictly positive continuous 
function on (0, oo). This means that the random variables ^'^(O) are independent and 
have Poisson distribution with rate p^{em). Similarly, the distribution of minus particles 
is (inhomogeneous) Poisson with density p_{em), where p_(x) is some strictly positive 
continuous function on (0, oo). 
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Notation 

Next we introduce main definitions and give their intuitive interpretation. Our inter- 
pretation concerns the situation of large density hmit when, instead of point particles on 
a lattice, there is a continuous media of infinitesimally small particles of two types, the 
particles move with fixed velocities — a+ < 0, a_ > and have initial densities p±{x) 
correspondingly. That is there are no fluctuations. Define the functions 

M^{r)= I p-{y)dy and M+(r) = / p+{y) dy for r > 0. (1) 



We interprete M±(r) as the cumulative mass of plus (minus) particles on the distance 
less than r from zero. Under above assumptions on p± we see that the functions M±(r) 
are strictly increasing on (0, +oo) and, therefore, the inverse functions r±{M), defined by 
the equation 

M±(r±(M)) = M 

exist and are strictly increasing. For example, the function r+ = r+(M) defines the 
interval (0,r+) where the mass of plus particles equals M. Then the function 

, , r_(M) + r+(M) , , 

T{M) := ^ !^ +y I ^2) 

defines the time interval (0,T(M)) during which mass M of plus and mass M of minus 
particles annihilate. The function T{M) is also strictly increasing on [0, +oo) and is 
invertible. Denote its inverse function by M{T). The place where the latter of these 
particles meet 

r+(M(T)) - a+T = -r_(M(T)) + a^T = (3{T) (3) 

is the coordinate of the boundary at time T. Excluding from the system the terms 
that are linear in T, we get 

^(T) = r^{M{T))^^ r_(M(T))- 



Scaling limit for the stochastic model 

Here we return to the stochastic particle model and formulate the main result. 

Theorem 2.1 For any fixed r > the following convergence in probability holds 

where the function (3 : R+ — > R zs deterministic and has the following explicit form 

-a+r_(M(r)) + a_r+(M(r)) 



Corollary 2.1 Consider the homogeneous case p^{y) = p-, y < 0, p+{y) = p+, y > 0. 
All functions defined above are linear: M_{r) = p^r, M+{r) = p+r, r±{M) = M/p±, 

r(M) = M ^ 

a_ + a+ 

M(T) = T _ + , 
P- +P+ 
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and, hence, the phase boundary Pslre moves with an asymptotically constant velocity: 

Pe{re )^/3{t) = t — — — 



r ■ 



P- +P4 

P+ + P- 



3 Proof 

Our plan is to show that the hmiting behavior of Ps{t) iii the stochastic model corresponds 
to the deterministic evolution described in (II])-(|3]). To do this we need some control over 
the random fluctuations in the limit e — i- 0. Now we fix M and consider the following 
random variables Ai± = Ai±{M), i = 0, 1,2, (we will prove that they are of the order 

1. Denote by Qe,± = Qe,±{M) the number of (±)-particles which were at time t = 
correspondingly in the intervals 

/; = (0,r+(M))nZ, and 1° = (-r_(M), 0) n Z,. (4) 

Define Ao,± = Qe,± - Ms^^. 

2. All particles among them will be annihilated during time t = T{M)e~^ except of 
the number Ai^± of them. 

3. Define A2^± as the number of plus and minus particles which were not at time t in 
the intervals (jl]) but were annihilated during time t = T{M)e~^. 

This control can be achieved by use of exponential bounds for some families of events. 
The proof uses some ideas from [1]. 

Definition 3.1 We say that a family of events A = {As}^^q has a property of exponential 
asymptotic sureness (e.a.s.) if there exist constants /C^ > 0, > 0, £^ > such that 
for all e < Eji, the following inequality holds 

P(A,)>l-/C^exp(-g^£-i). 

In the sequel, for breavity, we say sometimes that the event A^, has probability ex- 
ponentially close to one. We will use the following fact: if two sequences A = {As}^yQ 
and B = {i?e}^^Q have the property e.a.s., then this property holds also for the sequence 

C = {An 

It is helpful to enumerate the particles at time somehow with the only condition 
that 

■ ■ ■ < x^(0) < x^(0) < x^(0) < < x^(0) < a;^(0) < xj{0) < ■ ■■ 

Denote by g_(l) and the indices of plus and minus particles of the first annihilating 
pair. One can assume that if some plus (minus) particle jumps to a site where there are 
several minus (plus) particles then it annihilates with the minus (plus) particle having 
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minimal index. Let ai be the time moment when the first annihilation occurs. Since 
particles move independently, their order can change in time, so, in general, (i)(0) 7^ 
a;J~(0) and x'^^^^i^iS^) 7^ 2:^(0). Similarly, we define q-{m) and 5+(m) as the indices of 
the particles of the m-th annihilating pair and (Tm as the time moment of the m-th 
annihilation. 

Fix some M > 0. Let N, = [Me~^]. Consider the iV^-th pair of annihilating particles, 
•^q (ATg) •'^q+{N^)- '^^^ Hiain idea is to prove that for small e the random time a]\f^ 
is close to the value T{M)e~^ and the random coordinate x~ -^(O) G is close to 
— r_(M). In more precise terms, it is sufficient to prove that for any small fixed positive 
numbers xq, xi, ^_|_ with probability exponentially close to one (as s ^ 0) the following 
holds: 

(a) the moment a^^ belongs to the time interval {to{M,e),ti{M,e)), where 

to{M,e)^{T{M)-xo)e-\ 

ii(M,£) = (T(M) + xi)£-i; (5) 

(b) the starting point of the minus particle x~_^j^^-^{0) belongs to the set (— r_(M) — 
C-,-r_(M) + C-)nZ,; 

(c) similarly, the starting point of the plus particle a^^(^^)(0) belongs to the set (r+(M) — 

C+,MM) + c+)nz,. 

Let us prove the theorem assuming that the above statements (a)-(c) are proved. 
Recall that PeicNe + 0) = ^q_(N^){^Ne) — ^q+{N^)i'^Ne)- Individual motion of a minus 
particle is a simple random walk on Z^ with the mean drift «_£ = (A_ — //_) e, so 
applying the upper bound of the large deviation theory, we get that for any fixed i e N, 
s > and Sq > with probability exponentially close to one 

x~{se''^) — x^{0) e ((«_ - 5o) s, (q;_ + 60) s) . 

In fact, even stronger result holds: for fixed S2 > si > and do > the family of events 
{D^}, where 

-De = { x7(s£-i) - a;7(0) e ((«_ - Sq) s, («_ + So) s) , Vs e [si, S2] } , 
has a property of e.a.s. Together with (a) this gives 

^liN.)M - V(iv.)(0) e ((«- - So) {T{M) - xo) , {a. + So) {T{M) + xi)) 

with probability exponentially close to one. Combining the latter statement with the 
statement (b) we conclude that with probability exponentially close to one 

^q.iN,)M e (a-T(M) -r_(M) -7,a_T(M) -r_(M) +7) 

where 7 = "y{So, Xq, Xi, C_) > can be made arbitrary small, i.e., 

7(^0, Xo, xi, C_) ^0 as max((5o, xq, xi, (_) 0. 
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Using ([2]) we see that 



r_{M) + a_T{M) 



r_{M) + a 



r_{M) +r+(M) 



r_(M) 



+ r+(M) 



a;_ + «_(. 



and, hence, 



/3.(aivJ^-r-(M) 



a- + a+ 



+ r+(M) 



in probabihty as e — i- 0. 

To finish the proof of theorem we need only to check that 



This corresponds to continuity property of the border on the macroscopic time scale r = 
T. To establish this fact we should take into account that: 1) due to the drift assumption 
{a± > 0) the random sequence {<7m+i — cfm, rn G N} admits uniform exponential esti- 
mates for the tails of the distribution functions of o"m+i — cr^ (we refer the reader to [6] 
for the corresponding techniques); 2) in finite microtime t the displacements of walking 
particles have the order 0{e) while in finite macrotime r their displacements have the 
order 0(1). We omit the details. 

To prove the statements (a) we need the following main lemma. Denote by A/"- (0, tm{M, e)) 
a set of minus particles that collide with plus particles on the time interval (0, tm{M, e)). 

Lemma 3.1 For any sufficiently small X2, X3 > the following events 



have the probabilities exponentially close to one. 

Lemma 13.11 follows from Lemmas 13.21 and 13.31 Lemma 13.21 deals with the initial 
distribution of particles and Lemma 13.31 controls the deplacements of minus and plus 
particles. 

Lemma 3.2 Let t/i < 2/2 < and < zi < Z2. Then for any 6 > the following families 
of events have probabilities exponentially close to one: 



F, = {\Af_{OMM,e))\ < (M - xs)^-^} , 
= {|Ar_(0,ti(M,£))| > (M + X3)£-i} 




Re 



{the number of plus particles sitting at time t = in the set 
{zi, Z2) n is between [f^^ P+{y) dy — 6)e~^ 
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Lemma 3.3 For any 5i > each family of events 

Ae = {all particles xf{0) G (-r_(M) + 6i , r+(M) - 6i) n collide 
with particles of opposite sign till the time moment t{M)e^'^}, 
= {on the time interval t E (0, se^^) none of minus particles, started at 
t = from the set (— oo, — r_(M) — 6i) fl Z^, collides with any 
plus particle, started at t = from the set (r+(M) + 6i, +oo) fi Z^}, 

satisfies the e.a.s. property. Moreover, fix any y, x > and consider the following subsets 
ofZ,: 

S- = (-00, -y - x) , = {-y, 0) , = (0, y), S:^ = {y + x, +oo) . 
Define the events 

V^r = {on the time interval t G (0, se~^) none of minus particles, started at 

t = from S2 , will meet some minus particle, started from the set Si}; 

Us = {on the time interval t G (0, se~'^) none of plus particles, started at 

t = from S2 , will meet some minus particle, started from the set Si}. 

Then the families of events {V^} and {U^} have the e.a.s. property. 

Proofs of Lemmas 13.21 and 13.31 are based on standard probabilistic methods [B] and 
are omitted. Let us explain now how using these two lemmas one can get, for example, 
the upper bound for \J\f_{0,to{M,e))\ in Lemma [3A1 

Firstly, we include in this bound all minus particles starting at t = from the 
set (— r_ ( M(T(M) — xq) ) — ^5, 0) where ^5 > is small and will be fixed later. By 
Lemma [3.21 there is no more than (M„ (r_ ( M(T(M) — xq) ) + ^5) + Sq) of such par- 
ticles e.a.s. for small 6q > 0. 

We should add to this bound all minus particles that started at t = from the set 
(—00, — r_ ( M(T(M) — Xq) ) — ^5) and annihilated in the time interval (0, to{M, e)) with 
some plus particles. We will show now that with probability exponentially close to one the 
number N°{0,to{M,e)) of such minus particles can be estimated as C6~'^ where c > is 
any prefixed small constant. Indeed, by Lemma 13731 (again in the sense of e.a.s.) the minus 
particles in question can annihilate only by colliding with some plus particles, started at 
if: = from the set (0, r+ ( M(T(M) — xq) ) + 65). By Lemma [3.21 the number of the plus 
particles in this set is bounded by (M+ (r_|_ ( M(T(M) — xq) ) + ^5) + Sq) e~^. From this 
bound we should exclude plus particles which was annihilated in collisions with minus 
particles started from the set (— r_ (M(T(M) — xq) ) + ^5,0), since by the part "V^" of 
Lemma [3.31 during the time interval {0,to{M,e)) the latter minus particles will go ahead 
of the minus particles started from (—00, — r_ ( M(T(M) — xq) ) — ^5). By Lemma 13. 2[ 
initially there was no less than (M_ (r_ ( M(T(M) - xq) ) - ^5) - 6^) e'^ particles in the 
set 

(-r_(M(T(M)-xo)) + 55,0). 
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So using the mean value theorem from analysis we get 



e ■ iV°(0, to(M, e)) < (M+ (r+ ( M(T(M) - xq) ) + ^5) + ^e) 

- (M_ (r_ ( M(T(M) - xo) ) - ^5) - ^e) 
= M(T(M) - Xo) + M;(^i)55 + ^6 

- ((M(T(M) - Xo) ) - M:(^^2)55 - 5e) 
< S^iWM'Jc + WM'Jc) + 25^. 

Hence, in the sense of e.a.s., 

6-\Ar.{0,to{M,e))\ < {M_{r_{M{T{M) - Xo)) 

+ 6,) + 5e) + 55(||m:||c + ||m;||c;) + 2^6 

= M(T(M)-xo) + M:(^3)55 
+ S,i\\M'_\\c+\\M'Jc) + ^5e 

< M{T{M) - Xo) + 5,i2\\M'4c + ||M;||c) + 36^. 

It follows from ([2]) and assumptions on p± that M'{t) > k for some k > 0. Therefore, 
M{T[M) — Xo) < M — kxQ. Given xo > we are allowed to chose positive constants ^5 
and ^6 cis small as we like. So, finally, we get that with probability exponentially close to 
one the following estimate holds 

|Ar_(0,to(M,.))|< (m-^).-^ 

Lower bound for \J\f-{0,ti{M,e))\ can be obtained in a similar way. 

To get the proof of statement (b) one should combine Lemma 13.31 with the following 
lemma. 

Lemma 3.4 For any X5 > 0, the event 

H, = {x;iO) e (-(l + X5)r_(M),0) GAr_(0,ti(M,£))} 

has the property of e.a.s. 

Proof of this lemma uses arguments similar to the proof of Lemma [XT! The statement 
(c) is just a symmetric modification of the statement (b). 
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